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Abstract 

We compute the vector analyzing power (VAP) for the elastic scattering of trans- 
versely polarized electrons from protons at low energies using an effective theory of 
electrons, protons, and photons. We study all contributions through second order 
in E/M, where E and M are the electron energy and nucleon mass, respectively. 
The leading order VAP arises from the imaginary part of the interference of one- 
and two-photon exchange amplitudes. Sub-leading contributions are generated by 
the nucleon magnetic moment and charge radius as well as recoil corrections to the 
leading-order amplitude. Working to 0(E/M) 2 , we obtain a prediction for A n that 
is free of unknown parameters and that agrees with the recent measurement of the 
VAP in backward angle ep scattering. 
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1 Introduction 



The study of the vector analyzing power (VAP), A n , in polarized electron-proton 
scattering has recently become a topic of considerable interest in nuclear physics. The 
VAP is a time-reversal (T) odd, parity (P) even correlation between the electron spin and 
the independent momenta associated with the scattering process: 

A n ~ f r°Pp ii S 1 ,K a K , fi , (1) 

where S, P, and K (K') denote the electron spin, initial proton momentum, and incident 
(scattered) electron momentum, respectively. A non-zero VAP cannot arise at leading 
order in quantum electrodynamics (QED), but could be generated by new T-odd, P- 
even interactions involving electrons and quarks. Searches for such interactions have 
been carried out in neutron and nuclear /3-decay as well as nuclear 7-decays[TJ El E]. 
Indirect constraints may also be obtained from limits on the permanent electric dipole 
moments of neutral atoms under various assumptions regarding the pattern of symmetry- 
breaking jH El 13 Ej- The sensitivity of direct searches for T-odd, P-even interactions 
is generally limited by the presence of QED "final state interactions" (FSIs) that break 
the T-symmetry between initial and final states and give rise to non-vanishing T-odd, P- 
even observables. Uncertainties in theoretical calculations of these final state interactions 
would cloud the interpretation of a sufficiently precise T-odd, P-even measurement in 
terms of new interactions. Observations of T-odd, P-even correlations in nuclear 7-decays 
are consistent with theoretical calculations of QED final state interactions 0, while T-odd, 
P-even searches in neutron /3-decay have yet to reach the sensitivity needed to discern 
these effects. 

Recently, the SAMPLE collaboration has reported a non-zero measurement of the VAP 
in polarized, elastic electron-proton scattering [TO], making it the first non-zero result for 
any T-odd, P-even observable in any electron scattering process. The result has received 
widespread attention, as it differs substantially from the simplest theoretical estimate of 
QED final state contributions that neglects proton recoil and internal structure^J- While 
one might speculate that this difference reflects the presence of new physics, a more likely 
explanation lies in elements of nucleon structure omitted from the simplest treatments of 
QED FSIs. 

If so, then the SAMPLE result, as well as other VAP measurements that have been 
completed or are under consideration, could have important implications for the inter- 
pretation of other precision observables involving hadrons that require computation of 
QED corrections to the leading order amplitude. Such observables include the ratio 
of proton electromagnetic form factors obtained via Rosenbluth separation in elastic ep 
scattering[T2], higher-order "box graph" contributions to weak interaction observables [13 , 
or QED final state interactions in direct searches for T-odd, P-even effects. In each in- 
stance, a calculation of QED corrections requires a realistic and sufficiently precise treat- 
ment of hadronic intermediate states, particularly those arising in two-photon exchange 
amplitudes, .M 77 , or the analogous amplitudes involving the exchange of one heavy gauge 
boson and one photon. Since the leading QED contribution to A n arises from Im .M 77 , ex- 
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perimental measurements of the VAP provide an important test of theoretical calculations 
of A"l 77 needed for the interpretation of other measurements. 

At the same time, the VAP provides a new window on nucleon structure, as «M 77 
probes the doubly virtual Compton scattering (VVCS) scattering amplitude. In recent 
years, virtual Compton scattering (VCS) on the proton has become an important tool in 
probing the internal structure of the proton. VCS involves the coupling of one virtual 
and one real photon to a hadronic system. In the case of the proton, the VCS cross 
section is sensitive to the generalized polarizabilities of the proton, and its measurement 
should provide insight in the proton structure . In practice however, this cross section 
includes Bethe-Heitler (BH) amplitudes associated with radiation of a real photon from 
the electrons. Proper treatment of the cross section must therefore be taken in order to 
obtain a correct interpretation of the measurement. In contrast, the process involving 
the coupling of two virtual photons to the hadronic system is immune to background BH 
amplitudes and, thus, offers an alternative to VCS in probing the proton structure. 

With the aforementioned motivation in mind, we study the VAP in the framework 
of an effective theory of low-energy ep scattering. Since the SAMPLE measurement 
corresponds to kinematics close to the pion electroproduction threshold, we consider only 
the electron, photon, and nucleon as dynamical degrees of freedom. In this respect, our 
analysis corresponds to the use of heavy baryon chiral perturbation theory with the pions 
integrated out. To make the treatment systematic, we expand A n in powers of p/M, 
where p is either the incident electron energy (E) or mass (m) and M is the nucleon 
mass. Working to second order in p/M, we obtain all contributions to A n that arise 
uniquely from one-loop, two-photon exchange amplitudes and obtain a prediction that is 
free from any unknown parameters. We also write down the leading, non-renormalizable 
T-odd, P-even eepp operators whose intereference with A4 7 can generate a non-zero VAP 
and show that they contribute at O (p/M) 4 . 

We find that inclusion of all one- loop effects through 0(p/M) 2 in .M 77 as well as 
all terms in A^ 7 through this order is sufficient to resolve the disagreement between the 
SAMPLE result and the simplest potential scattering predictions. This resolution follows 
from several effects that occur beyond leading order in p/M: recoil corrections to the pure 
charge scattering result obtained in Ref. jllj . the nucleon isovector magnetic moment, 
and the proton charge radius. In the absence of dynamical pions, contributions from the 
nucleon polarizability arise at higher order than we consider here and appear unnecessary 
to account for the experimental result. Given that the incident electron energy E is of the 
same order as m^, we have no a priori reason to expect agreement of our computation with 
experiment. What it suggests, however, is that for this kinematic regime, pions play a less 
important role in the VVCS amplitude than one might naively expect. Future, low-energy 
A n measurements, taken over a broader range in q 2 and scattering angle than relevant to 
the SAMPLE measurement, would provide additional, useful tests of this conclusion. 

We also consider A n at forward scattering angles and energies somewhat higher than 
those of the SAMPLE experiment, since preliminary results for this kinematic domain 
have been reported by the A4 Collaboration at the MAMI facility in Mainz JH]- Although 
we would not expect our framework to be reliable in this kinematic regime, where the 
electron energy E is much closer to M, it is nonetheless instructive to compare with the 
Mainz preliminary results as a way of pointing to the physics that may be operative in 
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this domain. Indeed, we find substantial disagreement with the preliminary Mainz data. 
The culprit could be that going to the Mainz kinematics exceeds the limit of validity of 
our effective theory, that we must include additional dynamical degrees of freedom such 
as the 7r or A(1230) resonance, or both. Future studies using alternative methods such 
as dispersion relations may be needed to explore this kinematic domain. 

Finally, we also consider A n for polarized M0ller scattering. The VAP for this process 
has been measured by the E158 Collaboration at SLAC JH], and theoretical computations 
given in Refs. [T71 UHl HH] • Our computation agrees with these earlier A n (ee) calculations, 
providing a useful cross-check on our study of the VAP for ep scattering. 

Our discussion of these points is organized in the remainder of the paper as follows. In 
Section 2, we discuss general features of A n and our approach to the computation. Section 
3 provides details of the calculation. In Section 4, we give numerical results and discuss 
their significance, while Section 5 gives our conclusions. Technical details are provided in 
the Appendices. 



2 General Considerations 

We are interested in computing the VAP in elastic ep scattering: 

dcr T - da l 21m M* 11 M 1 



da^ + da l \M 7 \ 2 ' {) 



where dcru^ is the differential cross section for scattering of electrons with incident spin 
parallel (anti-parallel) to K x K' . In a phase convention where the single 7-exchange 
amplitude M. 1 is purely real, A n requires a non- vanishing imaginary part of yW 77 1 . To 
compute the latter, one must consider both the box and crossed-box diagrams of Fig. 1. 
Simple power-counting arguments indicate that the contribution to yW 77 arising from the 
leading-order ^p couplings is ultraviolet finite but infrared divergent. Thus, in general, one 
must also compute the contributions to A n arising from the bremsstrahlung diagrams of 
Fig. 2. As we show by explicit calculation in Appendix A, however, the bremsstrahlung 
contribution to A n vanishes identically, while IhlM 77 is infrared finite. The resulting, 
leading-order contribution to A n is 0(p/M) . 

Additional contributions to A1 77 arise from higher-order operators that couple one or 
more virtual photons to the proton and electron. We neglect the latter since they are 
suppressed by additional powers of the fine structure constant 2 . In contrast, the jp oper- 
ators are induced by strong interactions and have couplings of order e. In order to treat 
their contributions systematically, we adopt an effective theory framework since we can- 
not compute the operator coefficients from first principles in Quantum Chromodynamics. 
The natural framework for doing so is heavy baryon chiral perturbation theory (HB^PT), 
which provides a systematic expansion in powers of p/A x and p/M, where A x = AirF„ 
is the scale of chiral symmetry-breaking and p is an external momentum or mass with 

1 By lmA4 17 , we mean the coefficients of the various products of fermion bilinears, eTeNT'N, etc. 
that appear in the amplitude. 

2 For high energy scattering, these higher-order QED contributions may receive logarithmic 
enhancements ^5] . 
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Figure 1: (Color online) Two photon exchange diagrams. The wavy lines indicate virtual 
photons, while k (k f ) and p (p') denote the initial (final) electron and proton momenta, 
respectively. 



magnitude much less than M and A x . In the present case, where we integrate out the 
pions, we take p = E or m and use M as the heavy scale. For the kinematics of the 
SAMPLE experiment, E » m. Since there are no hard collinear infrared singularities in 
Im _M 77 , we may drop all power corrections involving the electron mass and obtain our 
result as an expansion in E/M. 

The leading terms in heavy baryon Lagrangian for nucleons and photons relevant to 
our computation are 

J_ = 
2 AT 



^N'y = B v iv ■ DB V + - ; - B v 



[vDY-D* 



B 



V 



+^e^F^v a B v S^B v - ^B v v^B v d x F» x + ■■■ (3) 

where B v is the field for a heavy proton of velocity v^, where = — ieA^, and 
where we have shown explicitly all jp interactions through 0(p 3 ). The latter arise from 
the subleading kinetic term in Eq. (J3J) as well as from the operators containing the 
field strength, F^ u . The coefficient \x = 2.793 is the proton magnetic moment, while C r 
determines the proton Sachs, or electric, radius: 

Cr = «L^ t = , (4) 

6 at 

where t = q 2 . The experimental value for (t 2 )e = 0.743 fm 2 [201 E] implies C r = 2.81. 
When included in the loop diagrams of Fig. 1, these interactions generate contributions 
to the ep amplitude M.^ and .M 77 through order (p/M) 2 relative to the leading term. 
To this order, operators associated with the nucleon polarizability (see Fig. 3e) do not 
contribute, as they occur at C(p 4 ) in C^-y when the pion is treated as heavy. 

Higher-order contributions to A n can also arise from effective T-odd, P-even eeNN 
interactions. The origin of such operators could be either physics that we have integrated 
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Figure 2: (Color online) Bremsstrahlung contributions. Labels are the same as in Fig. 1. 

out, such as contributions to .M 77 from ttN or A intermediate states , or explicit T-odd, 
P-even interactions arising from new physics. As shown in Appendix B, there exist no 
Hermitian, four-fermion operators at dimension six that contribute to A n . The lowest di- 
mension T-odd, P-even four fermion operators have dimension seven and would nominally 
contribute to A n at 0(p/M) 3 . We show, however, that contributions from these operators 
vanish to this order and first arise at 0(p/M) 4 . Since we truncate our analysis at two 
orders lower, we may neglect these operators and obtain a parameter-free prediction for 
the VAP. Nevertheless, we discuss these operators briefly in Section 4 when considering 
the possible size of neglected, higher-order contributions 3 . 

As we show in detail in Section 3, the leading one-loop contributions to A n - generated 
by two 0(p) 7p insertions in the VVCS amplitude (Fig. 3a) - are finite, non-analytic in 
p, and occur 0(p/M) , whereas those generated by the dimension seven T-odd, P-even 
operators arise at 0(p/M) A . Thus, the leading contributions are uniquely determined 
from the one- loop calculation. Similarly, contributions to M.^ involving one 0(p) and 
one 0(p 2 ) 7£> interaction (Fig. 3b, c) contribute to A n at 0(p/M), are also finite and non- 
analytic in p, and are unique to the loop calculation. The 0(p/M) 2 loop contributions 
arise either from two 0(p 2 ) ^p operators (e.g., two insertions of the nucleon magnetic 
moment operator, Fig. 3d) or one 0(p) and one 0(p 3 ) term (viz, the proton charge 
radius). We find, however, that the 0(p/M) 2 components of .M 77 arise only from the jp 
magnetic moment interaction as well as from recoil order terms in C^ 1 . Contributions 
to A4 77 from the proton charge radius vanish, though it does contribute to A n as a 
higher-order term in A / l 7 . 

3 For an earlier, phenomenological calculation that includes some of these higher order contributions, 
seeRef. [H]. 
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Figure 3: (Color online) Contributions to the VVCS amplitude appearing in Fig. 1. Open 
circles indicate the leading order jN couplings, while dark circles indicate higher-order 
couplings, such as the magnetic moment and charge radius. Shaded circle denotes that 
nucleon polarizability operator. 

3 Two-photon exchange 

The evaluation of four point functions for general kinematics does not readily lend itself 
to evaluation using standard Feynman parameterization in the loop integrals. Alternate 
methods for evaluating these integrals that do not rely explicitly on Feynman parameters 
have been worked out in Refs. j22l E3] arid have become standard. In the present case, 
where we are interested in backward angle scattering at nonzero q 2 , we would ideally like 
to use this formalism. However, the form of the heavy baryon propagator does not permit 
one to adopt the t'Hooft-Passarino-Veltmann formulation directly. 

We circumvent these difficulties by carrying out the computation with relativistic 
baryon propagators and expanding our result in powers of p/M. Doing so allows us to 
evaluate the loop integrals using the standard formulation of Refs. [211123] • It has been 
shown in other contexts [2E] that doing so allows one to recover the heavy baryon result 
so long as the external momenta are sufficiently small. Moreover, our loop results are 
entirely non-analytic in p and, thus, must match the corresponding non-analytic results 
obtained with heavy baryon propagators. To the order of our analysis, there exist no 
four fermion operators that could account for differences between relativistic and non- 
relativistic treatments of A n . 

The one-loop .M 77 is nominally infrared singular and must, therefore, be regulated 
with an IR regulator such as a photon mass. On general grounds, the regulator dependence 
should be cancelled by a corresponding dependence of the bremsstrahlung contribution 
to the spin-dependent cross section. As is well known, such a cancellation occurs for 
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unpolarized scattering cross section. In Appendix B, we work out the corresponding 
bremsstrahlung contribution to A n and show that it vanishes identically. Consequently, 
InxA/fyy must be IR regulator-independent. 

In general, the amplitude .M 77 depends on each of the eleven integrals obtained in 
Ref. [24J. The imaginary part, however, depends on only four: 

0,(1,2,3) = ^ln(A)e(s-(m + M) 2 ) 

VA sX 2 
C (l,3,4) = C (1,2,3) = C 



B (l,3) = 7r—Q(s-(m + M) 2 ) (5) 
s 

where the three labels associated with the B and Cq functions indicate which propagators 
are used for the two-point and three-point integral as discussed in Appendix C, A is the 
photon mass, and 

A = s 2 - 2s(M 2 + m 2 ) + (M 2 - m 2 ) 2 (6) 

These integrals have been previousely computed in Refs. |2"4*ll2l)] (In they are obtained 
by the use of dispersion techniques). The D Q and Co loop integrals diverge as A — > 0, but 
the combination 

2C + D t = \n(—)e(s - (m + Mf) (7) 
v A —st 

is finite in this limit and is the only combination of Dq and Cq integrals that is so. As 
such, the two-photon contribution to A n must only contain terms proportional to this 
combination or to the B integral. 

In evaluating the loop contributions to A n , it is most efficient to identify the terms 
in .M 77 that generate the correlation of Eq. (0) by carrying out the Dirac algebra in the 
interference term Im.M 77 A4* before evaluating the momentum integrals 4 . After carrying 
out the momentum integration, the contribution from the box diagram of Fig. la is 



/box/u* (47ra) 2 16m7T 2 (47ra) g 



4(M 2 - m 2 - 3s)M 2 R + k[(QR + 2)A - ((m 2 — M 2 — s)R + 2s)t] 
+ A(llm 4 - 2(13M 2 + 8s)m 2 + 15M 4 + lis 2 + UM 2 s)t 



+ 4s(2m 4 - (5M 2 + 4s)M 2 + 3M 4 + 2s 2 - 3M 2 s)t 2 ] 

A + ts , 9 
- {k 2 + Ak + 2)B 



(2C + D t) 



J This procedure introduces no ambiguities because Im Al 77 is finite to the order of our analysis. 
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s, £, and u are the Mandelstaam variables, k — pL — 1 is the nucleon anomalous magnetic 
moment and 

M 2 



i?- 1 = t 



AM 2 



To obtain the result consistent with our power counting, we expand Eq. 
of p/M up to second order relative to the leading term: 



(9) 
in powers 



ImM h ™M; 



+ 
e 



(4™) 2 327r 2 amM ^ 

t ^w^ 2 [{E 2 -m 2 + t/A) + §- I + ^] e F ^ K « K ? 



In 



A(E 2 



m 



2M 1 4M 2 
2\ Ii\j2 



2E/M + (2E 2 -m 2 )/M 



3E 2m 2 k 2 32{E 



m 



2\2 



+ t 2 /2 + 10(E 2 -m 2 )t 



Ak 
M 2 
(s- 



m 2 - E 2 
(m + M) 2 



M 2 

k 2 + Ak + 2 
M 2 



A(E 2 -m 2 )+t 



(E 2 -m 2 ) + 



(10) 



where the 0-function arises from the integrals 2Cq + 2Dot and Bq. Note that we have re- 
tained the m-dependence purely for illustrative purposes, as m << E for the experiments 
of interest here. The corresponding contribution from the crossed-box diagram can be 
obtained by crossing symmetry with the replacement s — > u. In this case, the 0-function 
vanishes, so only ImA^^° x A^* contributes. 

In the expression ()1U|). the terms that go as powers of E/M or m/M but do not 
contain factors of k or C r arise purely from recoil effects. The proton charge radius 
contributes solely via A4 T Although it also contributes to the absorptive part of A4 y7 , the 
resulting terms do not contribute to the spin-dependent correlation of Eq. Including 
the magnetic moment, charge radius, and recoil-order terms in M. 1 along with the loop 
contributions in Eq. (fTU|) leads to the following expression for the VAP: 



2atm 



VE 2 -m 2 [(E 2 -m 2 + t/A) + <§ + 
~A(E 2 - m 2 ~ 



S-K xK' 



In 



-t 
2m 2 



R+ M + M 2 



k 2 + Ak + 2 



2E/M+(2E 2 -m 2 )/M 2 

2 on(rp2 ^2\2 , -2 



k 32(E — m ) + t 12 + 1Q(E — m )t 



Ak 



M 2 



(E 2 



M 2 



m 



A(E 2 



+ t 



M 2 



E 



t 

+ 4 



X 



(8E 2 + 2t)R 2 + ^f + t 



AEt t + 2m 2 + 2k(* + 2m 2 ) + K 2 [t + A(m 2 - E 2 )}/2 



i -i 



M 2 



111 



Dropping all terms that go as powers of E/M, m/M, or t/M 2 yields the result obtained 
in Ref. JTTJ that was obtained for scattering from an infinitely heavy, point-like proton. 
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Figure 4: (Color online) VAP vs energy for fixed scattering angle, 9 = 146.1°. The dashed 
blue line is the leading order result, and the solid red line shows the full calculation. The 
SAMPLE result HH] is also shown at E = 192 MeV. 

4 Results and Discussion 

The expression for A n given in Eq. (jllj) provides a parameter-free prediction for low- 
energy electron scattering. In Fig. 4 and Fig. 5, we plot A n as a function of energy for 
fixed laboratory frame scattering angles 9 = 146.1° (Fig. 4) and 9 = 30° (Fig. 5), while 
in Fig. 6 we show the VAP for fixed energy E = 192 MeV while varying 9. In call cases, 
the leading order calculation is shown for comparison. In Fig. 6, the relative importance 
of the recoil, magnetic moment, and charge radius contributions are also indicated. 

The result obtained in the SAMPLE measurement is also shown. While the leading 
order calculation over estimates the magnitude of A n by a factor of roughly four, inclusion 
of the higher-order terms considered here produces agreement with the experimental value. 
Interestingly, there appears to be scant evidence that dynamical pions or the A play a 
significant role in A n for this kinematic region (E = 192 MeV), despite one's expectation 
that they might. 

At higher energies, our result for A n cannot be considered reliable, since the conver- 
gence of the effective theory expansion breaks down for E ~ M. The A4 collaboration 
at Mainz has measured A n at E = 570.3 MeV and E = 854.3 MeV and 25° < 9 < 35°. 
Preliminary results for the higher energy VAP have been reported in Ref. |15j . A compar- 
ison with our computation indicates that the preliminary experimental values for forward 
angle scattering and higher energies are substantially larger in magnitude than we are 
able to obtain via the low-energy expansion to 0(E/M) 2 . Presumably, a resummation of 
higher-order contributions in E/M using non-perturbative techniques, such as dispersion 
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Figure 5: (Color online) VAP vs energy for fixed scattering angle, 9 = 30°. The dashed 
blue line is the leading order result, and the solid red line shows the full calculation. 
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Figure 6: (Color online) VAP vs scattering angle for the SAMPLE kinematics (E = 192 
MeV). The dotted black line gives the leading order result, the blue dashed line adds the 
recoil corrections, the green dash-dotted line adds the magnetic corrections, and the solid 
red line shows the full calculation through 0(p/M) 2 . 
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Figure 7: (Color online) Possible contribution from the dimension seven, T-odd, P-even 
operator 7 e % to the backward angle VAP (9 = 146.1°). 

relations, would be required to compute reliably A n in this domain [l"2 7 l2~? t l2~£ | l2TJ | . We 
would also expect that inclusion of nucleon resonances 5 and pions as explicit degrees of 
freedom would be needed to account for the experimental results. 

One indication of the possible strength of these higher-order contributions may be 
given by considering the T-odd, P-even dimension seven operators. As shown in Appendix 
B, there exist two d = 7 operators that could, in principle, contribute. From an explicit 
calculation, we find that only one of the two - Ol% - leads to a non-vanishing A n . Here, 
it is useful to consider the form of this operator for relativistic proton fields, N: 

O 7 e a N = ^^ea^ l5 + ll) u eN l5lll N . (12) 
Re-writing this operator in terms of the heavy fields B v leads to 

dl a N = -2^j^ea^ 5 CD+^) u eB v S;B v , (13) 

where is the nucleon spin. The contribution from Ol% to the interference amplitude 

goes as e^ ual3 S^VyVaK'p and, thus, vanishes. On the other hand, using the 
relativistic form of the operator, 7 e %, leads to the correlation e^ val3 S ^P u P' a K'p that is 
non-vanishing for P ^ P' . The resulting contribution to the VAP is 

A (j) = aC 7a t 2 \K\\K'\sm6 

n ~ 4tt M 2 [8M 2 E 2 + 2(2E + M)tM + t 2 ] ' 1 ' 

5 For recent studies that pertain to such contributions, see Refs. 32, 33*1 134|. 
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Figure 8: (Color online) Possible contribution from 7 e ^ to the VAP at 9 = 30°, given 
constraints on the operator coefficient Cj a implied by the SAMPLE result. 

a result that is 0(p/M) A . In short, the only heavy baryon operators that can contribute 
involve either fields with two different velocities (viz, B v and B v >) whose contribution 
requires non-zero proton recoil, or dimension eight operators involving the B v fields only 
and carrying an additional p/M recoil suppression. 

The SAMPLE result for A n allows for a non-vanishing, but small coefficient for the 
leading, higher-order T-odd, P-even operator. Using the relativistic operator Olfj for 
illustration and including the loop contributions through 0(p/M) 2 leads to Cj a = 3.07 ± 
6.64. Naive dimensional analysis would have suggested a magnitude for Cj a or order 
unity, so the SAMPLE results do not appear to imply the presence of any un-natural 
hadronic scale physics. We may now use this range for Cr a to estimate the possible size 
of higher-order effects at other kinematics. The resulting band is shown in Fig. 7 for 
backward angles (9 = 146.1°) and in Fig. 8 for forward angles (9 = 30°). For the Mainz 
measurement at E = 570 MeV and 9 = 30°, we find -2.0 < A$ < 0.7 ppm, while 
A l ° op = —0.64 ppm. Thus, one might expect the impact of the physics we have integrated 
out to grow in importance relative to the loop effects considered here as the energy of 
the beam is increased, and it appears reasonable to expect a magnitude of a few ppm 
at the Mainz kinematics. We caution, however, that the precise value obtained in our 
calculation is unlikely to be correct in this energy regime, where the convergence of the 
E/M expansion is slow at best. 

As a final comparison, we also consider A n in fixed target, polarized M0ller scattering. 
The VAP for this process has been measured at SLAC by the E158 collaboration [TB], and 
one expects results to be forthcoming in the near future. Calculations of this quantity 
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Figure 9: (Color online) Diagrams contributing to the VAP for M0ller scattering. 



have been performed by several authors [T7J As a cross-check on our VAP for ep 

scattering, we carry out the analogous calculation here. It can be performed completely 
relativistically without performing an expansion in electron energy. However, since we are 
now dealing with identical particles in the final state we need to compute the interference 
between tree diagams in Figure 9b and the box diagrams of Figure 9a. For the SLAC 
measurement, one has E = 46 GeV. Performing the calculation in the center of mass 
frame we obtain: 



dcr^ da^ 



dtt dn 



o? m Am 2 

n a r- SHI U\ 1 

8 t 2 u 2 y/s V s 
3(s - Am 2 ) (t(u - s + 2m 2 ) In i 



-t 



s — Am 2 ' 



(15) 



-u(t - s + 2m 2 ) In ( 



—u 



s — Am 2 ' 



) - 2(t - u)tu 



da^ da^ 

~dn + ln 



a 



2st 2 u 2 



(t 2 + tu + it ) + Am 2 (m 2 -t -u)(t 2 - tu + u 2 ) 



Our results are in agreement with those of Refs. [T71 HH1 HH] 6 . The resulting asymme- 
try is ploted in Fig. 10, and agrees with the corresponding figure in Ref. ^Hj (note that 
in Ref. [IHj, the VAP is plotted vs. cos# rather than vs. 9 as we do here). 



6 In Ref. 0(a 2 ) contributions arising from initial and final state radiation effects were also 

computed. The corresponding contributions for the ep VAP are smaller than the hadronic uncertainties 
arising at Q(a), so we do not consider them 
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Figure 10: (Color online) The M0ller VAP vs CM scattering angle at the E158 kinematics. 

5 Conclusions 

In this study, we have computed the low-energy, backward angle VAP using an effec- 
tive theory involving electrons, photons, and protons and have obtained a parameter-free 
prediction through 0(p/M) 2 . The VAP to this order is determined entirely by the imag- 
inary part of the interference between the two-photon exchange, one-loop amplitude and 
the tree-level one-photon-exchange amplitude. In the limit that M — > oo, our result 
exactly reproduces the VAP obtained in Ref. ^T] for scattering from a structureless, 
infinitely heavy proton that over predicts the magnitude of A n at the kinematics of the 
SAMPLE experiment. We find that inclusion of all contributions through 0(p/M) 2 leads 
to agreement with experiment and leaves little room for important effects arising from 
dynamical pions or nucleon resonances at these energies. The leading counterterm contri- 
butions arise at 0(p/M) 4 and are consistent with zero. Thus, the SAMPLE measurement 
provides no evidence for unusual hadronic physics effects at these scales. The data also 
constrain the magnitude of the counterterm coefficients to be of natural size, and lead 
one to expect the VAP as measured by the A4 collaboration at Mainz to be at most of 
the order of a few ppm. Given the range of validity of our effective theory, however, we 
cannot produce a reliable prediction for VAP at the Mainz energies. 

In this context, the results of the SAMPLE measurement have notable consequences 
for studies of weak interaction processes. In the case of both neutron /9-decay and parity- 
violating ep scattering, theoretical consideration of final state QED corrections to the 
leading-order weak amplitudes is important for the interpretation of various measurements 
[T3j . To the extent that these measurements involve relatively low lepton energies, an 
analogous effective field theory computation of one-loop graphs involving the exchange 
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of one weak vector boson and one photon should be reliable at the ~ 20% level relative 
to the size of other 0(a) corrections. Future, more precise measurements of the VAP at 
low-energies and over a range of angles would provide important tests of this provisional 
assessment. 

One might also ask how competitive the SAMPLE measurement is with other direct 
searches for new T-odd, P-even interactions. As discussed in Refs.jHEj, direct searches 
are most relevant in symmetry-breaking scenarios wherein parity is broken at or above 
the scale for the breakdown of T. Existing direct searches imply that olt ~ few x 10~ 3 , 
where is the ratio of a typical T-odd, P-even nuclear matrix element to those of the 
residual strong interaction. When translated into bounds on generic, dimension seven 
operator coefficients CV [under the normalization of Eq. (j!2)) ]. one obtains \Cj\ < 2. 
The sensitivity of the SAMPLE measurement is comparable. Given that conventional, 
hadronic final state effects that have been integrated out in our computation naturally 
imply a value of Cj a with a magnitude of order unity, it appears unlikely that one will 
be able to circumvent the corresponding theoretical hadronic uncertainties as needed to 
make the VAP a direct probe of new physics. On the other hand, low-energy studies of A n 
could provide important information for the theoretical interpretation of other precision, 
electroweak observables. 
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Appendix A: Bremsstrahlung Computation 



Here, we show that the Bremsstrahlung amplitudes corresponding to Fig. 2 give a 
vanishing contribution to the VAP. The amplitudes are: 



M a 
M b 
M c 
M d 



q z (K — iy — m l 



—i 



2 

1 + 75^ 



«(iT)M(p / )(ie)7 /t M(p) 



-«(A-0(^e)7^(A-)n(pO(^e)7^^^^(^e)7 a 6 QM (p) 



q- " ' ' ' ' (p — I) 2 — M 2 

Here, 1^ is the radiated photon momentum. The square of the invariant amplitude 

2 



M 



B 



M a + --- + M L 



(16) 



(17) 
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depends on ten different products of leptonic and hadronic tensors. The leptonic tensors 
are: 



t ad _ j ac 



1 - ^ +m ^§^^ v+ ^§±i^^ 



r 6d r 6c 



t cd T cc 

(IV (IV 

Kt = 1%, (18) 
The corresponding hadronic tensors are: 

HZ = Tr[(j/ + M)'f(1> + M)'f] 

TlllV TTflu 

n ab — n aa 

HZ = Tr[(j/ + M)>f(j> + M)^ ( ^_~ M ^ 2 l»n 

HZ = Tr[(j/ + M)>f(j> + M)Y ^ 2 + _^ 2 1% 

HZ = HZ 

HZ = HZ 

n bd — n ad 

HZ = -[(^^(^^^^MK^^^/]^ 

HZ = Tr\(6' + M) 7 " ^ " ^ + M) -f{6 + M) 7 ^ ^ + M ) yw* 

We now need to compute: 



M B = 

pol 



t aa tjuu i t ab Tj^f i r ac rjni/ , r ad tjH-v i r 66 ttM 1 ' 

L fMuH aa + L^ v H ab + L H ac + L^ v H ad + L^ v H bb 
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i r be ttV-v _j_ t bd tj^ , r cc rr/i^ , j cd Tj^f , t dd Tj^f 

' lj f iu n bc + Ll t JLU n bd + L ^' fJ CC + ^flU U cd + ^liu n dd 



+ h.c. 



+ h.c. 



(20) 



where the sum is over all polarizations of the radiated photon. We are only interested in 
the terms proportional to e a p 1 sS a k l3 k'^p 5 . First we investigate the momentum integrals: 

11 11 



2n-B 



dH 



(p' + I) 2 — M 2 (p' + /) 2 — M 2 (p> + I) 2 - M 2 (p - Z) 2 - M 2 
11 11 

+ 7 ^ W~„ TTTT + 



(p-l) 2 - M 2 (p-l) 2 - M 2 (k' + I) 2 - m 2 (k' + I) 2 - m 2 
+ •••] (21) 
We can evaluate the generic two point integral as defined by: 



• 2 r>t 2 2 2\ 

27r ±j(p ; m 1} m 2 ) = /J, 



A—n 



d n q 



1 



X 



(22) 



_g 2 + m 2 — ie (q + p) 2 + m 2 — ie_ 

We are only interested in the imaginary part of B. We find that above the physical 
treshold s = —p 2 > (mi + m 2 ) 2 this integral develops an imaginary part Ref. [3~T] : 



T D / 2 2 2n \Av^™im|) 

LmB{p ; m l , m 2 ) = tt- B(s — (mi + m 2 ) J 



(23) 



Evaluating the B functions for the kinematics involved here we find that none of the 
integrals of Eqn. (|21|) develop an imaginary part. As such evaluating the traces and 
performing the integration we obtain a rezult of the form: 

M B = f l (m,M )S) t ) u) + f 2 (m,M )S) t ) u)ie a , 3lS S a k p k^p s + h.c. 

= 2f 1 {m,M,s,t,u) (24) 

Hence, we find no contribution to A n . 



Appendix B: Local Operators 

As discussed in the text, we are interested in computing the contribution to the VAP 
from local, four fermion eeNN operators. The lowest dimension operators of this form 
have dimension six. First, we show by explicit calculation that all d = 6 operators give 
vanishing contributions to A n . The most general form for the d = 6 operators are 

a 2 

0% = W2 e(Ci + C 2l ,)eN{C[ + C ! 2l ,)N 
a 2 

0% = j^eiCs + CwWeNiC^ + CM-f (25) 
a 2 

°eN = lH e{C 6 + C 1 ^)tT^eN{C' s + CMa iaf N (26) 
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where we have used relativistic nucleoli fields iV (the corresponding argument carries 
over straightforwardly in the heavy baryon formalism). To make the above hermitian 
we require all the constants C l eN to be real. We now compute the interference of the 
amplitudes associated with these operators and the tree amplitude .M 7 , retaining only 
the desired structure e al 3^sS a p l3 K 1 K' s . The corresponding leptonic and hadronic tensors 
are 



^6a 



6c 

Tlil 

n 66 
n 6c 



TrW + m)(d + C 275 )^r^^ 4 », h "] 



Tr[^' + m)(C 3 + C^)i 



2 

1 + 75^, 



M 6 M* +h.c. = 



Tr[r + m)(C 5 + O^K"^^ + m) 7 1 
7V[(^ + m)(C( + C^ 75 )(^ + m)Y] 
Tr[(fi' + m)(C 3 + C'^ 5 )r(P + m)r\ 
Tr[(t>' + m)(C 5 + C^s)^^ + m) 7 "] 
(47To;)a 2 



tM 2 



L^aH^a) + L^H^b) + L^H^qc) + , h.c. (27) 



Evaluating the traces and keeping only the terms of interest we obtain 

M e M; + h.c. = il6 {C X C' X M - C A C' A m)e aPl5 Sy p K^K' & + h.c. (28) 

Since all the C's are real we see there is no contribution from dimmension six terms. This 
results is as expected, as the operators C?6a-c are even under both T and P. 

Now consider d = 7 operators. As for the d — 6 operators, all contributions from 
T-even P-even d = 7 operators will vanish. We may, however, write down two Hermitian 
T-odd, P-even d — 7 operators: 



Ol a N = ^- 3 C 7a e^CD + H) u N l5lf xN 

a 2 4- -> 

Ol h N = j^C 7b e l5lfl eN l5 a^(D + D) V N 



(29) 
(30) 



As before we evaluate the interference of the above with M.~ r The corresponding leptonic 
and hadronic tensors are: 



L 7a 



L 7b 



m 7 m; + h. c . 



= i- 



iTrW + m^a^qj-^tf + m)Y] 

TrW + m ) 75 yl^% + m)Y] 

Tr[(fl + m)'y 5 'f(j> + m)'f] 
iTr\(i>' + m) lb a» a q a (i> + m)Y] 
. (4:7ia)a 2 



tM 3 



C7aL 7 aHfj,u(7a) + C 7b L^H^ 7b ) + h.C. (31) 
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Figure 11: (Color online) Momentum routing for the 77 box graph integrals. 

Evaluating the traces we note that only the L^H^^^ contributes: 

M 7 M, + h.c. = 16( y^ W (32) 

We are intrested in the contribution such a term gives to the VAP. Keep only the leading 
piece of the tree amplitude we get: 

A ( 7) = ®C7a t 2 \K\\K'\sm9 

4vr M 2 [8M 2 E 2 + 2(2E + M)tM + t 2 } 1 ' 

Appendix C: Loop Integrals 

Here, we provide additional details about the computation of As noted in the text, 

the contribution from the crossed-box diagram vanishes, so we consider only ImAi^Ai* . 
We may express the latter in terms of the leptonic and hadronic tensors: 



l z — mr 2 



2M v ' 1 ((I + K' + p') 2 - M 2 
x [i e {l + r{l + Kf) lv + ^{l + K) s \u{p) 

x u(p)[ie(l + r(K - K'f) la + ^(K - K'yW) 

p (H i — % — % — % 

M^° X M* = J j^yL^ a q - j^- 2 ^ - R ^ ^ - K , )2 H^ a (34) 

where 

r = R-l (35) 
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We define the loop integrals from above as follows 
m 2 Da = [ dH 



m 2 D n = / dH 



m 2 D a p = I dH 



iir 2 D a p 1 = I dH 



1 


(I 2 - 


m 2 )(l + K') 2 \(l + K> + p') 2 


-M 2 ](l + K) 2 


(I 2 - 


m 2 )(l + K') 2 \(l + K> + p') 2 

/ \ ' / L V ' JT / 


-M 2 )(l + K) 2 


(P- 


l a lpl^ 


-M 2 ](Z + if) 2 


(P- 


m 2 ){l + K') 2 [{1 + K' + p') 2 


-M 2 ](/ + iT) 2 


(P- 


m 2 ){l + K') 2 {{1 + K> + p') 2 


-M 2 ](l + K) 2 



™ 2 D^s = — ■ ■ (36) 



In order to evaluate these integrals, we follow the methods of Refs. j23j El], and our 
notation follows that of Ref. [23]. To this end, we need to compute the following three 
point functions 



m 2 C (l,2,3) = J dH 
m 2 C (l,2,4) = JdH- { 



1 



(P -m 2 )(l + K') 2 [(l + K' + p') 2 - M 2 } 
1 



(P - m 2 ){l + K') 2 {1 + K) 2 
m 2 C (l,3,4) = j dH ' 



(P - m 2 ) [(I + K' + p') 2 - M 2 } (I + K) 2 
m 2 C (2,3,4) = JdH {l - RI)2[{1 - — + ^ )2 - M2]{1 - K)2 



and two point functions 

m 2 B {l,2) = J dH 



(37) 



(P - m 2 )(l + K') 2 
m 2 fi (l,3) = JdH {l2 _ m2m + K/+pl)2 _ M2] 

= J dH (P-J K l + K) 2 

- 2B ^ = I dH (l + Kni + K) 2 

m2Bo(2 ' 3) = J dH (l + K<m-rK' + p') 2 -M 2 ] 
m2Bo(3 ' 4) = I dH [(l + K> + P r-M 2 )(l + K) 2 



(38) 



For all the B,C and D integrals above we are interested only in the imaginary part. The 
only two-, three- and four-point integrals with non-vanishing imaginary parts are: 
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Im D = 



- t H^)^=Q{s-{m + Mf) 



ImC (l,2,3) 
ImC (l,3,4) 

Im ^0(1,3) 



(m + M) 2 ) 



Im[Co(l,2,3)]=C 
7T— 9(s- (m + M) 2 ) 



(39) 



2\2 



In the above A is the photon mass and A = s 2 — 2s(M 2 + m 2 ) + (M 2 — m 

Although space considerations preclude a complete delineation of the calculation here, 
it is instructive to consider in more detail the evaluation of one of the four-point integrals 
required. Specifically, we consider 

D a = p°D u + P2*D 12 + p«D 13 (40) 

For the kinematics considered here the Passarino and Veltman momenta and masses are: 

Pi — K mi — m 
p 2 =p m 2 = 

P3 = —p' TO3 = M 

Pi = —K' m 4 = 

We then have for the Im 



where 



h = rn\ 



m 



p\ = 2m 2 



h = ml-ml+pl-pl = (s-M 2 

r 2 2 2,2 r 

J3 = m 2 -m A -p 4 +p 5 = -/ 2 , 

and where the inverse of the momentum matrix X is: 

-l 



m 2 ) 



X 



-l 



v\V2 pm 
P2P3 

pm P2P3 p\ 




A+ts 



3M 2 +m 2 -s-t 
A+ts 

-t 2(M 2 +s+t)m 2 -~(s+t-M 2 ) 2 - 
t(A+ts) 
M 2 -m 2 _ 1 
A+ts t 



M 2 -m 2 +s 
A+ts 
M 2 -m 2 1 
t 



A+ts 
s 

A+ts 



-1 J 



(41) 



where 


( Du 
Im D12 
V D 13 


j = X _1 Im 


f R20 \ 

R21 
{ R22 J 


(42) 


-R20 = 


^ [/iA> + 0,(1,3,4) - 


Co(2,3,4)] = 


: i(2ZW + C ) 




R21 = 


\ [/ 2 A> + C (1, 2,4) - 


C (l,3,4)] = 


:^[2D (s-M 2 -m 2 )-Co] 




R22 = 


^[/ 3 A> + C (1,2,3)- 


Co(l,2,4)] = 


-- l -[-2Do(s-M 2 -m 2 ) + C ] 


, (43) 



(44) 



(45) 
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After performing the necessary algebra we obtain 

£> [2((m - M) 2 - s)((m + Mf - s) + (m 2 — M 2 + s)t] - 2C (s + M 2 — m 2 ) 



lm[D n ] 
lm[D 12 ] 



2(A + ts) 

Doim* + (t- 2(M 2 + s)) + (M 2 - s)(M 2 - s - i))m 2 + 2C (m 2 - M 2 ) 



2(A + ts) 



Similar steps are required in evaluating the other four-point integrals. 
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